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We present a novel scheme for the appearance of Stochastic Resonance when the dynamics of a 
Brownian particle takes place in a confined medium. The presence of uneven boundaries, giving rise 
to an entropic contribution to the potential, may upon application of a periodic driving force result 
in an increase of the spectral amplification at an optimum value of the ambient noise level. This 
Entropic Stochastic Resonance (ESR), characteristic of small-scale systems, may constitute a useful 
mechanism for the manipulation and control of single-molecules and nano-devices. 
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Stochastic Resonance (SR) describes the counterintu- 
itive phenomenon where an appropriate dose of noise is 
not harmful for the detection or transduction of an in- 
coming, generally weak signal, but rather of constructive 
use in the sense that a weak signal becomes amplified 
upon harvesting the ambient noise in metastable, non- 
linear systems [l|. Since its first discovery in the early 
eighties SR has been observed in a great variety of sys- 
tems pertaining to different disciplines such as physics, 
chemistry^ engineering, biology and biomedical sciences 
0, B H S, S, U 0, S, E [l3|- The list of models and ap- 
plications is still growing. In particular, SR has found 
widespread interests and applications within biological 
physics. 

The research on SR has primarily been focused on sys- 
tems with purely energetic potentials. However, in sit- 
uations frequently found in soft condensed matter and 
biological systems, particles move in constrained regions 
such as small cavities, pores or channels whose pres- 
ence and shape play an important role for the SR- 
dynamics [lO||, sometimes even more important than the 
well-studied case of energetic barriers in such systems 
[111 . Il2l . Il3l . llj|. In this work, we demonstrate that ir- 
regularities in the form of confining, curved boundaries, 
being modeled via an entropic potential, can cause noise- 
assisted, resonant-like behaviors in the system under con- 
sideration. Confinement, an inherent property of small- 
scale systems, can thus constitute an important source 
of noise-induced resonant effects with interesting appli- 
cations in the design and control of these systems. 

The phenomenon of SR is rooted on a stochastic syn- 
chronization between noise-induced hopping events and 
the rhythm of the externally applied signal, that taken 
alone is not sufficient for the system to overcome a po- 
tential barrier. In the first place, noise enables system 
transitions and it is in fact responsible for the observed 
signal amplification and the emergence of certain degree 
of order. In the earliest and basic manifestation of SR, 
the synchronization of the random switches of a Brown- 
ian particle with a periodic driving force were observed 
for a bistable potential. Moreover, potentials of this type 
are not only found in systems with energy barriers, as 



FIG. 1: Schematic illustration of the two-dimensional struc- 
ture confining the motion of the Brownian particles. The 
symmetric structure is defined by a quartic double well func- 
tion, cf. Eq. ^, involving the geometrical parameters Lx, 
Ly and b. Brownian particles are driven by a sinusoidal force 
F{t) along the longitudinal direction and a constant force G 
in the transversal direction. 



they may also arise due to the infiuence of entropic con- 
straints. Particles diffusing freely in a confined medium 
such as the one depicted in Fig. [l] may give rise to an 
activation regime when a constant force G in the trans- 
verse direction is imposed. We will show that the com- 
bination of forcing and the presence of entropic effects 
deriving from the confinement and the irregularity of the 
boundaries give rise to an effective bistable potential that 
exhibits the signatures of Stochastic Resonance. 

The dynamics of a particle in a constrained geometry 
subjected to a sinusoidal oscillating force F{t) along the 
axis of the structure and to a constant force G in the 
transversal direction can be described by means of the 
Langevin equation, written in the overdamped limit as 



dr 



-Ge; - F[t)e, + yf/kBT ^(t) , 



(1) 



where r denotes the position of the particle, 7 is the 
friction coefficient, el and e^ the unit vectors along the 
X and y-directions, respectively, and £,{t) is a Gaussian 
white noise with zero mean which obeys the fiuctuation- 
dissipation relation {^i (t) ^j (t')) —2 Sij 5{t — t') for i,j — 
x,y. The explicit form of the longitudinal force is given 
by F{t) = i^osin(rit) where Fq is the amplitude and fi is 
the frequency of the sinusoidal driving. 



In the presence of confining boundaries, this equation 
has to be solved by imposing reflecting boundary condi- 
tions at the walls of the structure. For the 2D structure 
depicted in Fig. [H the walls are defined by 



WI{X) = Ly 



2 L,, 



L-. 



-Wuix), (2) 



where wi and Wu correspond to the lower and the up- 
per boundary functions, respectively, Lx corresponds 
to the distance between bottleneck position and posi- 
tion of maximal width, and Ly refers to the narrow- 
ing of the boundary functions and b to the remain- 
ing width at the bottleneck, cf. Fig. [ij Consequently, 
2w(x) = Wu{x) — wi{x) gives the local width of the struc- 
ture. This choice of the geometry is intended to resemble 
the archetype setup for SR in the context of a double 
well potential. In fact, in the limit of a sufficiently large 
transversal force, the particle is in practice restricted to 
explore the region very close to the lower boundary of the 
structure, recovering the effect of an energetic bistable 
potential. For sake of a dimensionless description, we 
henceforth scale all lengths by the characteristic length 
Lx, i.e. X = x/ Lx, y = y/Lx which implies b = b/Lx and 
wi — w\/Lx — —Wu, time by r = jLx /fceTR, the cor- 
responding characteristic diffusion time at an arbitrary, 
but irrelevant reference temperature Tr, i.e. t — t/r and 
fi ~ Ut, and force by i^R — jLx/t, transversal force 
G = G/Fji and a longitudinally acting, sinusoidal force 
F{i) — F{t)/Fji. In the following we shall omit the tilde 
symbols for better legibility. In dimensionless form the 
Langevin-equation JT]) and the boundary functions ([2]) 
read: 

dr 



-Ge-y - F{t)ex + Vi^^W , (3) 

w\{x) = —Wu{x) — ~w{x) — ex — 2ex — 6/2 , (4) 



dt "-^y 



where we defined the aspect ratio e = Ly/Lx and the 
dimensionless temperature D = T/Tr. 

In the absence of a time-dependent applied bias, i.e. 
F{t) = 0, it has been shown by a coarsening of the de- 
scription [ISl . Ilg Il7| that the Langevin equation ^ can 
be reduced to an effective ID-Fokker-Planck equation, 
reading in dimensionless form 



^^ll-S -■'<--'- 



where 



V{x,D) = -D In 



— smh 



Gw{x) 
D 



(5) 



(6) 



and the prime refers to the derivative with respect to 
X. This equation describes the motion of a Brownian 
particle in a bistable potential of entropic nature. Re- 
markably, the effective potential does not only depend 
on the energetic contribution of the force G, but also 



on the temperature and the geometry of the structure 
in a non-trivial way. Notably, for the vanishing width 
at the two opposite corners of the geometry in Fig. ^ 
this entropic potential approaches infinity, thus intrinsi- 
cally accounting for a natural refiecting boundary. It is 
important to emphasize that this bistable potential was 
not present in the 2D Langevin dynamics, but arises due 
to the entropic restrictions associated to the confinement. 
In general, after the coarse-graining the diffusion coeffi- 
cient will depend on the coordinate x as well, but since 
in our case |w'(a;)| "C 1, this correction can be safely 
neglected, cf. Ref. [H, [13, [H, [H, [l^l • 

It is interesting to analyze the two limiting situations 
that can be obtained upon varying the value of the ra- 
tio between the energy associated to the transversal force 
and the thermal energy. For the energy-dominated case, 
i.e. Gw{x)/D > 1, Eq. dH) yields V{x) = -Gw{x) (ne- 
glecting irrelevant constants). In the opposite Hmit, i.e. 
for Gw{x)/D <C 1, the corresponding entropic potential 
function reads V{x,D) = —D hi[2w{x)]. 

Two-state approximation.- It is instructive to analyze 
the occurrence of stochastic resonance in the context of 
the two-state approximation. For a potential V{x) with 
barrier height AV the escape rate of an overdamped 
Brownian particle from one well to the other in the pres- 
ence of thermal noise, and in the absence of a force, is 
given by the overdamped Kramers rate |2l|, [22, l23| , read- 
ing in dimensionless units. 



MD) ^ V^"(^n..)l^"(w)| ^^p 



-AV 

15~ 



(7) 



where V" is the second derivative of the effective po- 
tential function, and with a;inax and a;min indicating the 
position of the maximum and minimum of the potential, 
respectively. 

For the potential given by Eq. ^ and the shape de- 
fined by Eq. |[2|) , the corresponding Kramers rate for tran- 
sitions from one basin to the other reads, in dimensionless 
units. 



Ge v/2 sinh(G&/iJ) s{n\i[G{b + 2e) / D] 
'^^''^-- sinh^[G(5 + 2e)/i^] ■ (') 

Spectral amplification.- The occurrence of Stochastic 
Resonance can be detected in the spectral amplification 
77 |23l- It is defined by the ratio of the power stored in the 
response of the system at frequency VL and the power of 
the driving signal, and which for the periodically driven 
two-state model, cf. Ref. [l||, is given in dimensionless 
units as 



1 



V = 



44P) 



1)2 44(L>) + f72 



(9) 



Next we demonstrate the occurrence of the resonance 
in the spectral amplification that signals the phenomenon 
of ESR. We demonstrate that ESR is neither a pecuHarity 
of the two-state approximation nor of the equilibration 
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FIG. 2: (color online) In dimensionless units, the dependence 
of the spectral amplification rj on noise level D for different 
driving frequencies, for the transversal force G = 1.0, the 
driving amplitude -Fb = 10~ , and for the width function 
w{x) = — ea;* + 2ea;^ + 0.01 with the aspect ratio e = 1/4. The 
solid lines correspond to the ID modelling, cf. Eq.lfTO)l and 
Eq. lfl2|l . whereas the dashed lines correspond to the two-state 
approximation, cf. Eq.(l9]). 



assumption used to derive the effective ID-Fokker-Planck 
equation. 

In order to study the appearance of Stochastic Res- 
onance we analyzed the response of the system to the 
appUed sinusoidal signal in terms of the spectral ampli- 
fication T]. In the presence of an oscillating force F{t) in 
the x— direction there is an additional contribution to the 
effective potential in Eq. ^ and the ID kinetic equation 
in dimensionless units reads 



dt 



l[D'-£^iV'i^.,D)-Fit))P 



(10) 



The numerical integration of the ID kinetic equation 
(fT0| was done by spatial discretization, using a Cheby- 
shev collocation method, and employing the method of 
lines to reduce the kinetic equation to a system of or- 
dinary differential equations, which was solved using a 
backward differentiation formula method. This results in 
the time-dependent probability distribution P{x,t) and 
the time-dependent mean value, defined as 



WO)./xP(x,,)d.^ 



(11) 



In the long-time limit this mean value approaches the pe- 
riodicity of driving [22| with angular frequency fi. After 
a Fourier-expansion of {x{t)) one finds the amplitude Mi 
of the first harmonic of the output signal. Hence, the 
spectral amplification 77 for the fundamental oscillation 
reads: 



Ml 



F^ 



(12) 



The comparison of the ID modeling and the two-state 
approximation in terms of the spectral amplification rj, 



cf. Eqs. ^ and (fT2l) . demonstrates the capability of 
the two-state approximation for small driving frequencies 
and amplitudes, cf. Fig. [21 

ID modelling vs. precise numerics (2D).- In order to 
check the accuracy of the description we compared the 
results obtained by the ID modelling with the results of 
Brownian dynamic simulations, performed by integration 
of the overdamped Langevin equation JlJ , for a 2D struc- 
ture (see Fig. [J) whose shape is described by Eq. ^. In 
our case we have used the aspect ratio e = 1/4 and the 
bottleneck width b = 0.02. The simulations were carried 
out by the use of the standard stochastic Euler-algorithm. 

Fig. [3] depicts the dependence of the spectral ampli- 
fication r] on the noise strength for different values of 
the driving frequency, the driving amplitude and the 
value of G. It is important to point out that the re- 
sults obtained from the ID modelling using the ID- 
approximation (lines) are in excellent agreement with 
the numerical simulations of the full system (2D) (sym- 
bols) within a small relative error. This agreement is 
due to the fact that the considered potential function is 
a smooth function (|w'(a;)| ^ 1), and in this situation 
our employed ID-approximation is expected to become 
very accurate [15, 17, 18]. 

Fig. [3K shows the dependence of the spectral ampli- 
fication 77 on the noise strength D for various driving 
frequencies at a fixed transversal force and forcing am- 
plitude Fq. Here, we observe an increase in the spectral 
amplification which gives rise to the finding of the effect 
of Stochastic Resonance in the presence of entropic barri- 
ers. As observed for the usual "energy-dominated" SR []| 
the resonance peak is more pronounced as the applied an- 
gular frequency Jl of the input signal decreases. Similarly, 
Fig. \Sh depicts how ESR depends on the strength of the 
transversal force G. Interestingly, the maximal amplifi- 
cation increases upon decreasing the strength G of the 
transversal force. However, the presence of the transver- 
sal force G is crucial for observing a non-monotonic be- 
havior of the spectral amplification with increasing noise 
level D. In the limit of G ^ 0, the spectral amplifica- 
tion increases monotonically with decreasing noise level 
and tends asymptotically to l/il^ at small driving am- 
plitudes Fq. Finally, Fig. [S]: displays the dependence of 
the spectral amplification 77 on the noise strength D for 
various amplitudes of the driving force Fq at a fixed value 
of the transversal force and driving frequency. Both, the 
amplification of the signal and the optimal value of the 
noise at which it occurs seemingly increase as the driving 
amplitude decreases. 

This novel ESR effect is characterized by the appear- 
ance of a maximum in the spectral amplification as a 
function of the noise strength D, just as in conventional 
energy-dominated SR [l||. However, in biological and 
practical systems the temperature which controls the 
strength of the thermal noise is not a readily variable 
parameter. Our results suggest that, for a given tem- 
perature, a proper choice of the externally controlled pa- 
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rameters (i.e. the nature of the driving force, i.e., its 
amplitude and driving frequency and the strength of the 
transversal force) might bring the system into an opti- 
mal regime where confinement and noise mutually inter- 
play to boost noise-assisted transport inside a corrugated 
structure. 



In conclusion, we have elucidated a new mechanism 
leading to the appearance of noise-induced resonant ef- 
fects when a Brownian particle moves in a confined 
medium in the presence of periodic driving. The con- 
strained motion impedes the access of the particle to cer- 
tain regions of space and can be described in terms of 
a bistable entropic potential. The activated dynamics of 
the particle in this effective potential then results in a co- 
operative effect between noise and external modulation, 
yielding an Entropic Stochastic Resonance. The effect 
detected is genuine for small-scale systems in which shape 
and fiuctuations are unavoidable factors ruling their evo- 
lution. The advantageous possibilities of ESR on what 
concerns optimization and control may provide new per- 
spectives in the understanding of systems at the scales 
of micrometers and nanometers and open new avenues in 
their manipulation and control. 



FIG. 3: (color online) In dimensionless units, the dependence 
of the spectral amplification r] on noise level D for various 
values of the quantities, the driving amplitude Fo, driving 
frequency £7, and the transversal force G. The symbols cor- 
respond to the results of the Langevin simulations for the 
two-dimensional structure with the shape defined by the di- 



mensionless function w{x) 



-ex* + 2ex^ +0.01 with the 



aspect ratio e = 1/4, whereas the lines are the results of the 
numerical integration of the ID kinetic equation l|10p. 
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